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Silicon nanotube is constructed by rolling up a silicene, i.e., a monolayer of silicon atoms forming a two- 
dimensional honeycomb lattice. It is a semiconductor or an insulator owing to relatively large spin-orbit interac- 
tions induced by its buckled structure. The key observation is that this buckled structure allows us to control the 
band structure by applying electric field E z . When E z is larger than a certain critical value E a , by analyzing 
the band structure and also on the basis of the effective Dirac theory, we demonstate the emergence of four 
helical zero-energy modes propagating along nanotube. Accordingly, a silicon nanotube contains three regions, 
namely, a topological insulator, a band insulator and a metallic region separating these two types of insulators. 
The wave function of each zero mode is localized within the metallic region, which may be used as a quantum 
wire to transport spin currents in future spintronics. We present an analytic expression of the wave function for 
each helical zero mode. These results are applicable also to germanium nanotube. 
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INTRODUCTION 

Carbon nanotube is one of the most fascinating materi- 
als. There is a variety of nanotubes from metal to insu- 
lator depending on how it is constructed by rolling up a 
graphenelll \M- Similarly, silicon nanotube may be con- 
structed by rolling up a silicene 13j-|6|], a monolayer of silicon 
atoms forming a two-dimensional honeycomb lattice. Sili- 
con nanotubes have already been manufactured[7-9]. Almost 
every striking property of carbon nanotube is expected to be 
transferred to this innovative material since carbon and silicon 
belong to the same family in the periodic table. Nevertheless 
there exists a major difference. A large ionic radius of silicon 
induces a buckled structure|5[, which results in a relatively 
large spin-orbit (SO) gap of 1.55meV. Accordingly, silicon 
nanotube has a finite gap and it is always a semiconductor or 
an insulator. One might think it is an ordinary bund insulator 
just as most carbon nanotubes are. See FigQ]for an illustration 
of carbon nanotube and silicon nanotube. 

In this paper we reveal an amazing property of silicon nan- 
otube thanks to this buckled structure when it is placed in ex- 
ternal electric field. Analyzing the band structure in the pres- 
ence of electric field E z perpendicular to the nanotube axis, 
we demonstrate that silicon nanotube is actually a topological 
insulator. When E z is beyond a certain critical field E a , we 
find four zero-energy modes to emerge in the bulk band gap 
and form metallic regions along a nanotube. A silicon nan- 
otube is made of three different phases, the topological insu- 
lator region, the bulk insulator region and the metallic region 
separating them. 

Topological insulator OlOJ 1 1111 is a new state of quantum 
matter characterized by a full insulating gap in the bulk and 
gapless edges topologically protected. These states are made 
possible due to the combination of the SO interaction and the 
time-reversal symmetry. The two-dimensional topological in- 
sulator is a quantum spin Hall (QSH) insulator with helical 
gapless edge modes lll2ll . which is a close cousin of the integer 
quantum Hall stat e. Q SH insulator was proposed by Kane and 
Mele in graphene lll 311 . However, since the SO gap is rather 
weak in graphene, the QSH effect can occur in graphene only 



(a) Carbon nanotube (b) Silicon nanotube 




FIG. 1: (a) Carbon nanotube and (b) silicon nanotube. The lattice 
is distorted due to a large ionic radius of a silicon atom and forms a 
buckled structure in silicon nanotube. 



at unrealistically low temperature lll4l 11511 . Our finding is that 
it is materialized naturally in silicon nanotube. 

As we have stated, there emerge four zero-energy modes in 
silicon nanotube under uniform electric field. They are helical 
edge modes of a topological insulator and propagate along a 
nanotube: They transport only spins without charges. These 
observations are supported by the effective Dirac theory. We 
construct explicitly the wave functions describing the four he- 
lical zero modes. In conclusion, we are able to realize a dis- 
sipationless spin current along a silicon nanotube by applying 
uniform electric field. It may be used as a quantum wire in 
future spintronics. 



SILICENE AND TIGHT-BINDING MODEL 

The band structure of a silicon nanotube is obtained sim- 
ply by imposing a certain periodic boundary condition to a 
silicene sheet provided the diameter is large enough. Silicene 
consists of a honeycomb lattice of silicon atoms with two sub- 
lattices made of A sites and B sites. The states near the Fermi 
energy are it orbitals residing near the K and K' points at op- 
posite corners of the hexagonal Brillouin zone. We refer to 
the K or K' point also as the K^ point with the valley in- 
dex r\ = ±1. Due to the buckled structure the two sublat- 
tice planes are separated by a distance, which we denote by 2£ 
with £ — 0.23A. It generates a staggered sublattice potential 
oc 2£E z (x, y) between silicon atoms at A sites and B sites in 
external electric field E z (x, y). 

The silicene system is described by the four-band second- 
nearest-neighbor tight binding model lll6ll . 
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where c ia creates an electron with spin polarization a at 
site i, and (i,j) / ((i,j)) run over all the nearest/next-nearest 
neighbor hopping sites. The first term represents the usual 
nearest-neighbor hopping with the transfer energy t — 1.6eV. 
The second term represents the effective SO coupling with 
Aso = 3.9meV, where a — (<J x ,<7 y ,a z ) is the Pauli ma- 
trix of spin, Vij = (dj x dj) / |d{ x dj| with dj and dj the 
two bonds connecting the next-nearest neighbors. The third 
term represents the Rashba SO coupling with Ar = 0.7meV, 
where fa = ±1 for the A (B) site, and d^ = dy/ |dy |. The 
forth term is the staggered sublattice potential term. The same 
Hamiltonian as can be used to describe germanene, which is a 
honeycomb structure of germanium |pl 1 1 611 . where various pa- 
rameters are t — 1.3eV, Aso = 43meV, Ar = 10.7meV and 
£ = 0.33A. 

By diagonalizing the Hamiltonian ([T]) under uniform elec- 
tric field E z , the band gap A (E z ) of silicene is determined to 
be 



A (E z ) = 2 \£E Z - V s z X 



so 



(2) 



at the K^ point, where s z = ±1 is the electron spin. The gap 
A (E z ) closes at E z — ?]s z E c with 



E ct = X S o/£ = 17meV/A, 



(3) 



where it is a semimetal due to gapless modes. It has been 
shown[6] that silicene is a topological insulator for \E Z \ < 
E'er, while it is a balk insulator for \E Z \ > E cr , as illustrated 



in FigfJ] Hence a topological phase transition occurs between 
a topological insulator and a band insulator as E z changes. 

The topological insulator is characterized by one of the 
following two defining properties! lOJ ll 111 . (1) The topolog- 
ical insulator has a nontrivial topological number, the Z2 
index II 1311 . which is defined only for a gapped state. (2) There 
emerge gapless modes in the edges [Fig|2jb)]. These two 
properties are closely related one to another. Indeed, the rea- 
son why gapless modes appear in the edge of a topological 
insulator is understood as follows. When a topological insu- 
lator has an edge beyond which the region has the trivial Z2 
index, the band must close and yield gapless modes in the in- 
terface. Otherwise the Z2 index cannot change its value across 
the interface. 

The above criteria cannot be applicable to a nanotube as 
they are since it is intrinsically one-dimensional (ID). We 
overcome the problem by applying electric field and by creat- 
ing a domain with ID edges in the surface of a nanotube: See 
Fig|5]we derive later. 



SILICON NANOTUBE 

A silicon nanotube can be constructed by rolling up a sil- 
icene precisely just as a carbon nanotube is constructed by 
rolling up a graphene|[l] |2|]. There are a variety of ways how 
to rolling up a silicene as in the case of carbon nanotubes. It 
is specified by the chiral vector, 



L = niai + U2SL2, 



(4) 



where ai and a2 are the basis vectors of the honeycomb lattice 
with integers ri\ and ri2. Let us take the coordinate (x,y) with 
the x-axis parallel to the chiral vector L (the circumference 
direction) and the y-axis orthogonal to it (the nanotube axis 
direction). The nanotube's circumference is given by 



|L| = a\/rf + n\ — mn2, 



(5) 



where a = 3. 86 A is the lattice constant. Let us denote the 
conjugate momemtum as (hk x , hk y ). We impose the periodic 
boundary condition ip(x+L, y ) = ip(x,y). It makes the wave 
vector k x discrete, i.e., 

k x = 2nj/L, i = 0,l,---,2n-l, (6) 

where 2n is the number of silicon atoms per unit cell. The 
wave vecto 
louin zone, 



wave vector k y remains continuous within the ID first Bril- 



ir/T < k y < ir/T, 



(7) 



where T is the length of the primitive translation vector in the 
y direction. The ID energy bands are given by the straight 
lines specified by © and (O in the (k x ,k y ) space. The sim- 
plest nanotube is the armchair type, where L — ^/3na. 



(a) Topological insulator (nanoribbon) [£ z |<£cr 




(a) Topological insulator (nanotube) \E Z \ <E C r 



(b) Band insulator (nanoribbon) |£z|>£cr 
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FIG. 2: (Color online) (a) The band structure of silicene nanoribbon 
for \E Z \ < E a . There appear two bands crossing the gap since there 
are two edges, indicating it is topological insulator. When the Rashba 
interaction is neglected (Ar = 0), each band contains two-fold de- 
generate zero-energy states corresponding to up and down spins, (b) 
The band structure of silicene nanoribbon for \E Z \ > E a . All states 
are gapped, and it is bulk insulator. The horizontal axis is the mo- 
mentum, and the vertical axis is the energy in unit of the transfer 
energy t. 




(b) Coexistence (nanotube) 
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FIG. 3: (Color online) (a) The band structure of silicon nanotube for 
\E Z \ < E a . All states are gapped, and it is topological insulator, (b) 
The band structure of silicon nanotube for \E Z \ > E a . There appear 
four bands crossing the gap. They are the zero-energy states separat- 
ing the topological and band insulating states in the nanotube surface: 
See Fig(5] When the Rashba interaction is neglected (Ar = 0), each 
band contains two-fold degenerate zero-energy states corresponding 
to up and down spins. The horizontal axis is the momentum, and the 
vertical axis is the energy in unit of the transfer energy i. 



The energy spectrum of a silicon nanotube is simply given 



by 



~? V*V) — ^ (."'xj Ky)% 



(8) 



in terms of the energy spectrum £(k x , k y ) of silicene. The 
energy spectrum of a carbon nanotube can be gapless since 
graphene has a gapless band. There is a variety of carbon nan- 
otubes from metal to insulator depending on the way how to 
choose the chiral vector L. On the contrary, a silicon nan- 
otube has always a finite gap since silicene has a finite gap. 
The band gap is rather insensitive to the choice of L. 

We have calculated the band structure of an armchair sili- 
con nanotube based on the tight-binding model (Q} by chang- 
ing the uniform external electric field E z . The band structure 
is strikingly different for \E Z \ < E cr and \E Z \ > E cr , as illus- 
trated in Figf3] All states are gapped for \E Z \ < E cr . On the 
other hand, for \E Z \ > E CI , there appear four bands crossing 
the gap at the K and K' points. When the Rashba interaction is 
neglected (Ar = 0), each band contains two-fold degenerate 
zero-energy states corresponding to up and down spins. Even 



when Ar 7^ 0, the degeneracy remains unsolved at the K (K') 
point, though it is slightly resolved away from these points. 
The feature of the zero-energy states is highly contrasted with 
that in a silicene nanoribbon, where they emerge at k = n as 
in FigUta). 

We show the probability density of these four zero modes 
in Fig|U to which we have assigned the spin s z and the valley 
K M where they appear. We are able to make the assignment 
and determine the direction of the current based on the effec- 
tive Dirac theory: See a sentence below (flTt and Figj5] Here 
we explain them by using the gap formula (O and the band 
structure in Figj3jb). 

It follows from the formula (f2) that the electrons with spin 
s z = ±1 are gapless at the K± point for E z > and that 
the electrons with spin s z = =Fl are gapless at the K± point 
for E z < 0. Thus the peaks indexed by Kf appear in the 
upper half and those indexed by K| appear in the lower half 
of the nanotube. Recall that E z is the electric field between 
the A and B sublattices, which is opposite between the upper 
and lower halves of the nanotube. Furthermore, the direc- 
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FIG. 4: (Color online) The probability density of the zero modes 
in silicon nanotube. There appear four peaks at 6 — Oj, with 
sin 9j = ±E a /E, as indicates the emergence of four metallic states 
separating topological and band insulators: See Figj5] The holizontal 
axis is the circumference (0 < x < L) or the angle (0 < 6 < 27r). 



tion of the current can be determined by examining the dis- 
persion relation of the zero-energy mode in the band structure 
[FigHt))]. We find the up-spin electrons and the down-spin 
electrons propagate into the opposite directions in each peak. 
Namely the current is helical, that is, it is a spin current, as 
illustrated in Fig|5] It is a characteristic feature that helical 
zero modes appear along the edge of a topological insulator. 
We conclude that the zero-energy metallic state separates the 
topological and band insulating states in a silicon nanotube. 



DIRAC THEORY 

In order to explore deeper physics of the helical zero 
modes, we analyze the low-energy effective Hamiltonian de- 
rived from the tight binding model (H). It is described by the 
Dirac theory around the K n point as[ 16] 



with 



H„ = hv F (k x T x - i]k y T y ) + r/T z h n + £E z t z , (9) 



tin = -Asoc 2 - a\ R {k y (T x - k x (T y ) , (10) 



where t q is the Pauli matrix of the sublattice pseudospin, v F = 
^ at = 5.5 x 10 5 m/s is the Fermi velocity. 

We have shown numerically that there emerge four zero- 
energy states inside the bulk band gap as in Fig|3jb). The 
low-energy Dirac theory allows us to investigate analytically 
the properties of these helical zero modes. In so doing we 
set Ar = to simplify calculations. This approximation is 
justified by the following reasons. First all all, we have nu- 
merically checked that the band structure is rather insensitive 
to Ar based on the tight-binding Hamiltonian (Q}. Second, Ar 
appears only in the combination (k x ± ik y )\R in the Hamilto- 
nian (O, which vanishes exactly at the K± points. Third, the 
critical electric field E CT is independent of Ar as in ©. 
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FIG. 5: (Color online) An illustration of silicon nanotube under elec- 
tric field E > E a . There appear two topological insulator regions 
and two band insulator regions. They are separated by metallic states 
made of helical zero modes. A spin current flows in each metal- 
lic region as indicated. For instance, up-spin (down-spin) electrons 
propagate into the left (right) direction at 8 — 9\ . 



We take the y-axis parallel to the nanotube axis and the x- 
axis along the circumference. We may set k y =constant due to 
the translational invariance along the y axis. The momentum 
k y is a good quantum number. Setting 



&{x,y) = e lk y y <$>(x), 



(11) 



we seek the zero-energy solution, where Vt 7 (x, y) is a four- 
component amplitude. The particle-hole symmetry guaran- 
tees the existence of zero-energy solutions satisfying the rela- 
tion <pB (x) = i£,4>A (x) with £ = ±1. Here, <j>a is a two- 
component amplitude with the up spin and the down spin, 
(j>A = (<f>A~, 0a)* • Then the eigenvalue problem yields 



(12) 



H v 4>a{x) = E v ^4>a{x), 

together with a linear dispersion relation 

E ni = Vtyvpky. (13) 

The equation of motion for the component <j> 8 ^ (x) reads 

(Zhopd* + r}Sz\so - iE z (x)) 4> s l (x) = 0. (14) 

We apply the uniform electric field E perpendicular to the 
nanotube axis. 

The effective field for electrons in a silicon nanotube is 
given by 



E z {x) =Esm9 



(15) 



with 9 = 2nx/L. We take E > E CI with ©. We solve the 
equation 



E z {6)=ris z E a , 
Let us chose one solution and denote it as 

6 CI = Arcsin(E cr /E). 



(16) 



(17) 



We obtain two solutions 6\ = 9 CI and 62 = tt — 9 CI for r\s z = 
+1, and two solutions #3 = 7r + 9 CI and 64 = 2tt — 9 cr for 



rjs z = — 1, It implies that the zero modes at 8 = 9\ contains 
up-spin electrons (s z = +1) from the K valley (r/ = +1) and 
down-spin electron (s z = — 1) from the K' valley (jj = — 1), 
and so on, as illustrated in Fig]4] As we shall soon see, the 
sign £ is fixed at each Oj, as implies that up-spin and down 
spin electrons propagate into the opposite directions according 
to the dispersion relation ([T3V The current is helical, whose 
direction is determined by the sign of £. 
The equation of motion (TT~4T > is rewritten as 



—£hvpde — 2lEcos — — 

Lj Z 



' '#f(0) = o. 

(18) 



To construct the solution near Oj, we approximate 



cos- 



We can explicitly solve this as 



<^(0) = Cexp 



* 



2LIE 
Trhvp 



cos Oj cos ■ 



(19) 



(20) 



where C is the normalization constant. It satisfies the periodic 
boundary condition 4>a{9 + it) = 4>a{9) or <f>A (x + L) = 
4>a{x). The sign £ is to be chosen so as to make the peak of 
the wave function to appear at 9 — Oj, and we have 



(0) =C*exp 



2LIE . n . 0-0j 

— COS Oj cos — — 

nhv F ' J ' 2 



(21) 



We obtain £ = -1 at 9 = 0\ and = 9 4 , while £ = +1 
at 9 = 0-2 and 9 = 9%. As we have noticed in the sentence 
below ( fTTI i, the sign of £ determines the direction of the helical 
current at each Oj based on the dispersion relation (fT~3T >, which 
we have illustrated in FigfS] Furthermore, we have checked 
that the probability density \<ffi (9) | 2 agrees excellently with 
the result obtained based on the tight-binding model in Fig|4] 



CONCLUSION 

We have uncovered a salient fact of silicon nanotube that 
it is a topological insulator by studying its band structure in 



uniform electric field E z , When E z > E a , there emerges 
four helical zero modes propagating along the nanotube. They 
form metallic regions separating topological and band insu- 
lators. We have constructed the wave function of each zero 
mode based on the effective Dirac theory. Silicon nanotube 
may be an ideal material to transport spin currents. It may be 
used as a quantum wire in future spintronics. Our results are 
applicable also to germanium nanotube. 
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